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ABSTRACT 

The Sombrero galaxy, NGC 4594, contains the most numerous globular cluster system 
of any nearby spiral. It is an ideal candidate in which to study the globular clusters and 
contrast them with those in Local Group spirals. Here we present B and I imaging from 
the CTIO Schmidt telescope which gives a £eld-of-view of 3l' x 31. Using DAOPHOT 
we have detected over 400 globular clusters and derived their magnitudes, B-I colors 
and photometric metallicities. We have attempted to separate our sample into disk 
and bulge/halo globular cluster populations, based on location in the galaxy. There is 
some evidence that the disk population is more metal-rich than the bulge/halo globular 
clusters, however contamination, dust reddening and small number statistics makes this 
result very tentative. We hnd that the median metallicity of the bulge/halo globular 
clusters is [Fe/H] = -0.8. This metallicity is consistent with previous estimates based 
on smaller samples. It is also similar to the metallicity predicted by the globular cluster 
metallicity - galaxy luminosity relation. As with our Galaxy, there is no radial metallicity 
gradient in the halo globular clusters. This suggests that the spheriodal component of 
NGG 4594 did not form by a dissipational process. 


^ Based on observations obtained at Cerro Tololo Inter-American Observatory. CTIO is operated by AURA, Inc. under contract to the 
National Science Foundation. 

^Lick Bulletin 1350 
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1. Introduction 

Measurement of the abundances of Milky Way 
globular clusters (GCs) have provided important clues 
on the formation and chemical enrichment history of 
the Galaxy. Gombined with kinematics, a clear di¬ 
chotomy of metal-rich disk GGs and metal-poor halo 
GCs have been identified (e.g. Zinn 1985). The disk 
and halo GCs have a mean metallicity of [Fe/H] = - 
0.6 and -1.6 respectively (e.g. Ashman & Bird 1993). 
Globular clusters in our sister galaxy, M31, have been 
relatively well-studied so that spectroscopic metallic¬ 
ity measurements are available for about half of them 
(Huchra, Brodie & Kent 1991). A bimodal metallicity 
distribution is not as obvious for M31 but the analysis 
of Ashman & Bird (1993) indicates peaks at [Fe/H] 
= -0.6 and -1.5, i.e. similar to the Milky Way. The 
other Local Group spiral, M33, has about 30 known 
GCs, of which spectroscopic metallicities have been 
obtained for 22 (Brodie & Huchra 1991). The mean 
metallicity of [Fe/H] = -1.55 ± 0.37 is similar to that 
of Milky Way and M31 halo GCs. The small number 
of GCs makes it difficult to convincingly identify a 
population of disk GCs. 

Beyond the Local Group the nearest spirals are 
M81 (m-M = 27.8) and NGC 4594 (m-M = 29.73). 
Perelmuter & Racine (1995) estimated from imag¬ 
ing of M81 that it contains ~ 200 GCs. Spectra 
for 30 GCs have been obtained (Brodie & Huchra 
1991; Perelmuter, Brodie & Huchra 1995) which give 
a mean [Fe/H] = -1.48 ± 0.19. The Sombrero galaxy, 
NGC 4594, is a giant Sa galaxy (My = -22.0) with 
almost 2,000 GCs (Bridges & Hanes 1992). This is 
significantly more than either the Milky Way with 
Ngc = 160 ± 20 or M31 with Ngc = 350 ± 100 
(Harris 1996). It is also nearly edge-on {i = 84°) so 
that GCs projected close to the galaxy major axis can 
be used to define a plausible sample of disk GCs. Sep¬ 
arating disk and halo GCs in M31 (z = 38°) has been 
a major source of uncertainty. Thus NGC 4594 is an 
ideal candidate in which to study the GC system of 
another spiral galaxy. 

The distance to NGC 4594, from the average of the 
surface brightness fluctuation and planetary nebulae 
luminosity function methods, is 8.8 ± 0.4 Mpc (Cia- 
rdullo, Jacoby & Tonry 1993). At this distance GCs 
range in magnitude from B ~ 18-27, so that only the 
brightest few GCs are accessible spectroscopically by 
today’s large telescopes. A valiant effort was made 
recently by Bridges et al. (1996). After 7 hours of 


integration time over 4 nights of observing with the 
4.2m William Herschel Telescope, they obtained spec¬ 
tra of 34 confirmed GCs. Unfortunately their spectra 
were not of sufficient quality to determine a metal¬ 
licity for individual GCs. Summed together they de¬ 
rived a mean of [Fe/H] = -0.7 ± 0.3. An alternative 
approach, is to obtain GC colors from imaging stud¬ 
ies. Imaging has the advantage of being very efficient 
and photometric metallicities can be derived from the 
Galactic GC color-metallicity relation (Couture et al. 
1990). 

The GC system of NGC 4594 has been studied pho¬ 
tometrically by Harris et al. (1984) and Bridges & 
Hanes (1992). The first of these was a wide field-of- 
view photographic study in U and V bands. They 
examined the spatial distribution and estimated the 
GC specific frequency. The second study targeted 
three areas of the galaxy and imaged these in B and 
V bands with a 2.2 x 3.6 CCD. From the B-V colors 
of 131 GCs they derived a mean [Fe/H] = -0.81 ± 0.25 
for the GC system. Here we present 31 x 31 field- 
of-view CCD imaging of NGC 4594 in B and I bands. 
We measure B-I colors which are about twice as sen¬ 
sitive to metallicity as B-V colors (e.g. Geisler, Lee & 
Kim 1996). Our wide field-of-view imaging allows us 
to define spatially samples of ‘disk’ and ‘bulge/halo’ 
GCs, and examine their metallicity distributions. 

2. The Data and Globular Cluster Selection 

Broad-band B and I images of NGC 4594 were 
taken with the Cerro Tololo Interamerican Observa¬ 
tory (CTIO) Schmidt telescope. We used a Thomson 
1024 X 1024 array with a pixel scale of 1.84 /pixel. 
The images were obtained in 1994 February 12 and 
13. Reduction was carried out in the standard way 
(i.e. bias and dark subtraction, flat-fielding and sky 
subtraction). The total exposure times were 7200 and 
2400 secs for the B and I images respectively. After 
combining, the B image was calibrated using aperture 
photometry from the UBV catalog of of Longo & de 
Vaucouleurs (1983) which gave an rms accuracy of ± 
0.06 mag. Burkhead (1986) found that the bulge had 
a near constant color of B-I = 2.7 ±0.1 and this was 
used to calibrate the I image. Finally, the magnitudes 
are corrected for Galactic extinction of Ab = 0.12 and 
A/ = 0.05 (Bender et al. 1992). Typical photometric 
errors are 0.08 at B = 20.5 and 0.18 at B = 21.5. Our 
B-I colors are transformed into [Fe/H] values using 
the color-metallicity relation of Couture et al. (1990) 


2 



which is based on Milky Way GCs. This relation is 
calibrated from about 1/100 solar metallicity up to 
solar. The typical random error in this transforma¬ 
tion is ± 0.15 dex in [Fe/H]. From the final images, 
we subtracted a smooth elliptical model for the bulge 
using IRAF software. 

Contamination from foreground stars and back¬ 
ground galaxies is a concern for most GC imaging 
studies. Although not ideal, one method is to obtain 
similar exposure time images in some nearby ‘blank’ 
field and process these images in the same way as the 
actual data. This provides a statistical contamination 
correction. We do not have such a ‘blank’ field, so we 
attempt to estimate the contamination using star and 
galaxy surface densities from the literature. Before 
applying this method we will restrict candidate GCs 
to magnitudes of 17.5 < B < 22. The faint limit is a 
little brighter than the expected turnover of the GC 
luminosity function (i.e. B ~ 22.5) and was chosen 
to ensure that a color bias was not introduced. The 
bright limit is at least 3a from the turnover magni¬ 
tude assuming a dispersion characteristic of the Milky 
Way and M31 (Seeker & Harris 1993), while excluding 
the most obvious foreground stars. 

The number of foreground stars can be estimated 
from the Bahcall & Soneria (1981) model of the Milky 
Way. For our 31 x 31 image in the direction of NGC 
4594, their model predicts about 750 stars with 17.5 

< B < 22. Using figure 1 of Koo & Kron (1992), we 
estimate our image contains about 220 background 
galaxies down to B = 22. Thus the total contamina¬ 
tion in our image from stars and galaxies, before any 
color selection, is about 970 for magnitiudes 17.5 < 
B < 22. To compare this number with the number 
of objects in the B image, we have used DAOPHOT 
with a fairly conservative detection threshold of 5cr 
per pixel. This gives ~ 1400 objects with 17.5 < B 

< 22. Thus we estimate our image contains ~ 430 
bona fide GCs. Based on our detection limit and an 
expected turnover magnitude of B ~ 22.5, our sam¬ 
ple represents about 1/3 of the total GC population. 
So the total number of GCs is ~ 1300, which is at 
the lower limit of the range estimated by Bridges & 
Hanes (1992) i.e. 1900 ± 600. We note that although 
we have only detected a fraction of the total GC pop¬ 
ulation, there is no known correlation of GC metallic¬ 
ity with luminosity or mass (Huchra, Brodie & Kent 
1991) which might otherwise have bias our results. 

Although stars and galaxies cover a similar range 
of B-I colors as GCs, we know from other studies 


that the vast majority of GCs have metallicities -2 < 
[Fe/H] < -kO.5 (i.e. 0.15 < B-I < 2.33). So exclud¬ 
ing objects beyond these values will preferentially re¬ 
move stars and galaxies from the object lists. After 
applying this metallicity (color) selection and making 
restrictions in the DAOPHOT sharpness and round¬ 
ness parameters, our sample consists of 457 candidate 
GCs. This is similar to the ~ 430 objects left after a 
statistical correction for stars and galaxies has been 
applied. From these arguments we expect that the 
majority of objects in our final list are indeed GCs. 

Next we attempt to divide the sample into ‘disk’ 
and non-disk (which we will call ‘bulge/halo’) GCs. 
In our Galaxy (Gilmore & Reid 1983) and NGC 4594 
(Burkhead 1986) the scale height of disk stars is Z ~ 
1 kpc. Disk GCs in our Galaxy can reach distances 
above the plane of Z < 4 kpc. We have chosen to de¬ 
fine disk GCs in NGC 4594 as those with scale heights 
|Z| <4 kpc. The colors of the disk objects may be 
reddened by the galaxy’s dust lane. To minimize this 
effect we have excluded objects that are close to the 
central dust i.e., objects with major axis distances 
of less than 5 kpc of the galaxy center. Conversely 
bulge/halo GCs are simply defined as those with jZj 
> 4 kpc; it is not clear whether such GCs are associ¬ 
ated with the dominant bulge or the halo. Dividing 
the sample in this way we have 79 disk objects and 
378 halo objects. Due to projection effects some bu¬ 
gle/halo GCs will be present in the disk sample. From 
the surface density of halo GCs we estimate that be¬ 
tween 1 /3 and 1 /2 of the disk GCs are not associated 
with the disk. 

3. Results and Discussion 

In this section, we first discuss the GC sample of 
Bridges et al. (1996). Second, we examine the GC 
metallicity distribution in NGC 4594 and finally we 
search for spatial abundance gradients in this distri¬ 
bution. 

After superposing the position of the Bridges et al. 
34 spectroscopically-confirmed GCs on our CCD im¬ 
ages we managed to measure B and I magnitudes for 
24 (the others are presumably fainter than our lim¬ 
iting magnitude). The mean color is B-I = 1.9 ± 
0.2 (error on the mean), which corresponds to [Fe/H] 
= -0.7 ± 0.5. This is reassuringly consistent with 
the spectroscopic mean metallicity, for all 34 GCs, of 
[Fe/H] = -0.7 ± 0.3. This suggests that the transfor¬ 
mation from color to metallicity is reasonable. 
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In Fig. 1 we show the metallicity distribution for 
disk (|Z| < 4 kpc) and bulge/halo (|Z| > 4 kpc) ob¬ 
jects. The median metallicity of the bulge sample is 
[Fe/H] = -0.8. This is similar to that obtained by 
Bridges et al. (1996) and suggests that our sample 
is dominated by bona fide GCs. As pointed out by 
Bridges et al. the GC metallicity is more comparable 
to that found in giant ellipticals than other spirals. 
However this might be expected given the claimed 
correlation of GC metallicity with parent galaxy lumi¬ 
nosity (e.g. Brodie & Huchra 1991). Using the recent 
best-fit of Forbes et al. (1996) and the bulge luminos¬ 
ity of NGC 4594, the relation predicts a GC system 
mean metallicity of [Fe/H] = -0.6 ± 0.3 (rms error). 
The disk sample has a median metallicity of [Fe/H] 
= -0.5. As the disk sample will include a significant 
fraction of bulge GCs, the average metallicity of the 
disk sample is likely to be somewhat higher than - 
0.5. If we separate the 24 Bridges et al. (1996) GCs, 
with available B-I colors, into bulge and disk sys¬ 
tems, then their metallicities are [Fe/H] = -1.2 and 
-0.1 respectively. Thus there is some tentative evi¬ 
dence, from both our sample and that of Bridges et 
al. that NGC 4594 contains a population of disk GCs 
that are slightly more metal-rich than non-disk GCs. 
Although it is difficult to rule out preferential redden¬ 
ing of the GCs by dust, we do not think this has had 
a strong effect on the disk GC colors. First, we have 
excluded the inner galaxy regions with the obvious 
dust lane. Second, we do not find a strong radial GC 
color gradient in either scale height (see Fig. 2) or 
along the major axis as might be expected from cen¬ 
trally concentrated dust. The possibility that inner 
GCs are instrisically blue (and so removing a redden¬ 
ing trend due to dust) is very unlikely given studies of 
other galaxies (e.g. Harris 1991; Geisler et al. 1996). 

Abundance gradients, or lack of, provide useful 
constraints to galaxy formation models. For example, 
the absence of an abundance gradient in Milky Way 
halo GCs supports the view of Searle & Zinn (1978) 
in which the halo is made up of small protogalactic 
fragments, rather than arising from a pure dissipa- 
tional collapse. In the M31 GCs the situation is more 
uncertain. Huchra, Brodie & Kent (1991) found that 
there exists a upper envelope in GC metallicity with 
galactocentric radius but also pointed out that selec¬ 
tion effects and contamination by disk GCs made evi¬ 
dence for an abundance gradient inconclusive. In Fig. 
2 we show the spatial variation of GC metallicity for 
the disk and bulge/halo samples. In the lower panel 


of Fig. 2 we show GC metallicity versus projected 
galactocentric radius for the bulge/halo objects. As 
with the Milky Way GC system, there is no statisti¬ 
cally significant metallicity gradient. This constrains 
the role of dissipation in forming the spheroidal com¬ 
ponent of NGC 4594. In a dissipational collapse one 
expects a metallicity gradient in the direction of the 
collapse e.g., perpendicular to a disk. Accretion of 
material may reduce the strength of such gradients. 
Our Galaxy shows a weak trend of GC metallicity 
with scale height, where the slope is -0.084 ± 0.031 
dex kpc“^ (Armandroff 1993). In the upper panel of 
Fig. 2 we show the metallicity of disk GCs versus 
scale height from the major axis. The data show a 
hint of a weak trend with scale height, but it is not 
formally significant. 

4. Conclusions 

After separating a sample of ~ 450 globular clus¬ 
ters in NGC 4594 with photometric metallicities into 
spatially-defined ‘disk’ and ‘bulge/halo’ subsystems, 
we find tentative evidence for a difference in the 
metallicity of the two populations. The bulge/halo 
globular clusters have a median metallicity of [Fe/H] 
= -0.8 whereas the disk globulars have [Fe/H] > -0.5. 
Contamination, dust reddening and a small number of 
disk globular clusters make it difficult to confidently 
ascribe a higher average metallicity to the disk pop¬ 
ulation. The metallicity of the bulge/halo system is 
greater than that of the Milky Way ([Fe/H] = -1.6) 
or M31 ([Fe/H] = -1.5) halo globulars, however it 
is consistent with that expected from the correlation 
of globular cluster metallicity with galaxy luminos¬ 
ity. Like the Milky Way, we find no evidence for a 
metallicity gradient in the bulge/halo globular clus¬ 
ter population. This argues against formation of the 
spheroidal component by dissipational collapse. 
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Fig. 1.— Metallicity distribution of globular cluster candidates in NGC 4594. The upper panel (the ‘disk’ sample) 
shows objects with scale height |Z| < 4 kpc (excluding the inner galaxy regions). The median metallicity is [Fe/H] 
= -0.5 although there is clearly a large range in metallicities. The lower panel (the ‘bulge/halo’ sample) shows 
objects beyond 4 kpc from the galaxy disk. The median metalliticity is [Fe/H] = -0.8. There is tentative evidence 
that the disk globular clusters are more metal-rich than the bulge/halo globular clusters. 
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Fig. 2.— Spatial variation of globular cluster metallicity in NGC 4594. The upper panel shows the variation of 
metallicity with scale height for the disk sample (i.e |Z| < 4 kpc). There is no statistically significant gradient. The 
lower panel shows the radial variation of metallicity for the bulge/halo sample (i.e |Z| > 4 kpc). As with the Milky 
Way, there is no statistically significant radial metallicity gradient. 
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